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Abstract 

The  effect  of  the  choice  of  control  allocation  algorithm 
used  in  conjunction  with  a  nonlinear  adaptive  control 
law  is  examined.  In  particular,  an  existing  backstepping 
control  law  design  is  modified  to  allow  the 
incorporation  of  various  control  allocation  algorithms. 
Simulation  results  are  presented  for  two  separate 
adaptive  control  law  designs  and  for  a  non-adaptive 
two-loop  dynamic  inversion  controller.  With  the  control 
law  designs  fixed,  tracking  maneuvers  are  simulated 
using  four. separate  control  allocation  routines  -  Direct 
Allocation,  Discrete  Time  Direct  Allocation,  Pseudo- 
Inverse,  and  Weighted  Pseudo-Inverse.  System 
performance  is  then  examined  as  a  function  of  control 
allocation  method.  As  expected,  even  though  the 
remainder  of  the  control  law  design  remains  fixed, 
system  performance  is  directly  impacted  by  the  choice 
of  control  allocation  algorithm.  Further,  some  of  the 
important  control  law/control  allocation  interactions  are 
identified  and  their  effects  on  overall  system 
performance  analyzed. 

Introduction 

Basic  control  law  algorithms  and  control  allocation 
approaches  are  fi*equently  treated  as  separate 
disciplines.  Papers  on  control  laws  often  use  fairly 
simple  control  allocation  approaches  or  assume  the 
same  number  of  controls  as  controlled  variables. 
Papers  on  control  allocation  often  look  primarily  at  the 
ability  to  attain  commanded  moments.  Yet,  as  both 
flight  control  laws  and  control  allocation  techniques 
become  increasingly  complex,  the  possibilities  of 
adverse  or  unexpected  interactions  increase.  While 
there  have  been  some  efforts  that  have  looked  at  the 
advanced  control  law/control  allocation  problem  (e.g., 
References  1-5),  many  questions  still  exist.  This  is 
particularly  true  with  regard  to  potential  interactions 
with  adaptive  or  other  highly  nonlinear  control  laws. 

'Flight  Control  Engineer,  Member  AIAA 
Tlight  Control  Engineer,  Senior  Member  AIAA 
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Further,  for  some  advanced  control  algorithms,  there  is 
the  possibility  of  improvements  through  greater 
integration  between  the  control  law  and  control 
allocation  algorthms.  For  example,  in  the  case  of  an 
adaptive  control  law,  should  the  control  allocation 
parameters  remain  fixed  or  should  the  control  allocation 
be  updated  in  some  way  to  deal  with  errors  in  control 
effectiveness  caused  by  model  uncertainty  or  aircraft 
damage?  For  an  adaptive  controller  that  explicitly 
identifies  control  effectiveness,  there  may  be  a 
reasonably  straightforward  way  of  performing  the 
update  (although  there  are  a  number  of  practical 
difficulties  involved  such  as  effector  collinearities).  For 
a  direct  adaptive  controller,  the  situation  is  much  less 
clear.  Some  direct  adaptive  approaches  have  no  known 
way  of  extracting  control  effectiveness  changes.  Some 
other  approaches,  like  backstepping  adaptive  control, 
have  pseudo-estimates  of  control  parameters,  but  the 
parameter  estimation  is  done  to  meet  a  total  system 
stability  criterion  rather  than  to  optimize  the  error  in 
estimation.  As  a  result,  the  estimates  are  not 
guaranteed  to  converge  to  the  actual  values  over  time 
and  it  is  not  clear  whether  the  use  of  such  values  in 
control  allocation  would  be  beneficial.  This  is  further 
complicated  by  the  fact  that  different  inputs,  different 
flight  conditions,  or  different  levels  of  model 
uncertainty  may  have  substantial  effects  on  the  values 
of  the  estimates  and  the  resulting  control  law 
interactions. 

This  paper  examines  the  interactions  of  two  adaptive 
backstepping  control  law  designs  and  a  non-adaptive 
two-loop  Dynamic  Inversion  control  law  with  four 
different  control  allocation  approaches.  Dynamic 
Inversion  is  an  approach  based  on  canceling  the 
aircraft’s  natural  dynamics  and  substituting  desired 
dynamic  behavior^.  Backstepping^  is  a  Lyapunov 
approach  that  can  be  used  to  design  stable  adaptive 
control  laws  for  a  wide  range  of  systems.  The  basic 
concept  behind  backstepping  is  to  use  some  states  as 
virtual  controls  to  control  other  states.  The  primary 
benefits  of  this  type  of  controller  are  that  it  allows  a 
wide  array  of  nonlinearities  to  be  incorporated  in  the 
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controller  design,  and  it  has  proven  nominal  stability 
and  convergence  of  error.  Previous  work  in  the 
application  of  adaptive  backstepping  designs  to  the 
aircraft  control  problem^’^  have  assumed  the  same 
number  of  control  as  controlled  variables.  Concurrent  to 
the  work  presented  herein,  Reference  9  examines  the 
use  of  backstepping  control  with  redundant  control 
effectors.  While  Reference  9  uses  only  a  Weighted 
Pseudo-Inverse  solution  to  the  control  allocation 
problem,  this  paper  endeavors  to  begin  an  examination 
of  the  dynamic  interactions  between  various  control 
allocation  schemes  and  nonlinear  adaptive  controllers 
of  the  type  examined  in  References  7-9. 

If  there  are  no  limits  on  the  control  effectors,  then  the 
particular  control  allocation  algorithm  employed  is 
usually  considered  to  be  of  little  consequence.  That  is  to 
say,  any  legitimate  control  allocation  algorithm  will  be 
able  to  realize  the  desired  control  effect.  However,  if  an 
estimate  of  the  control  effectiveness  is  used,  different 
control  allocation  algorithms  could  indeed  produce 
different  responses  even  when  no  constraints  are  placed 
on  the  control  effectors.  Furthermore,  when  control 
effector  rate  and/or  position  limits  are  considered,  the 
control  allocation  algorithm  becomes  extremely 
important  in  determining  system  response  even  with 
perfect  knowledge  of  the  control  effectiveness 
parameters. 

There  have  been  many  recent  advances  in  control 
allocation  (e.g..  References  3,  5,  10,  and  11).  For  the 
initial  examination  of  the  effect  of  control  allocation  on 
system  behavior  under  the  control  of  a  nonlinear 
adaptive  control  law,  four  control  allocation  techniques 
will  be  implemented:  Direct  Allocation  Discrete 
Time  Direct  Allocation*  ^  Pseudo-  Inverse,  and 
Weighted  Pseudo-Inverse.  While  a  more  exhaustive 
array  of  control  allocation  algorithms  could  be 
implemented,  those  chosen  should  prove  sufficient  to 
demonstrate  some  of  the  important  control  law/control 
allocation  interactions  and  the  profound  impact  of  these 
interactions  on  system  performance. 

Aircraft  Model 

The  chosen  nonlinear  aircraft  dynamics  model  (Eq.  1) 
is  similar  to  that  used  in  References  7-9,  Here,  Aa  =  a  - 
a©  and  the  control  positions  are  measured  from  their 
trim  values.  The  seven  independent  control  effectors  are 
taken  to  be  the  left  and  right  horizontal  tails,  left  and 
right  ailerons,  collective  leading  edge  flaps,  collective 
trailing  edge  flaps,  and  collective  rudders.  Some  of  the 
key  simplifications  made  in  this  model  are  constant 
velocity,  no  lift  and  drag  effects  of  the  control  surfaces, 
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and  no  higher  frequency  dynamics  (e.g.,  actuators).  The 
control  surface  movements  are  however  restricted  by 
the  rate  and  position  limits  given  in  Table  1.  The 
nominal  model  parameters  are  given  in  Table  2  for 
trimmed  flight  at  30,000  feet  Mach  0.7.  Table  3  gives 
the  parameter  values  for  trimmed  flight  at  40,000  feet 
Mach  0.6.  Trim  values  for  the  horizontal  tails,  leading 
edge  flaps,  and  trailing  edge  flaps  are  -0.68°,  4.41°,  and 
5.56°  for  Table  1,  and  -0.97°,  -0.70°,  and  10.73°  for 
Table  2  -  those  not  specified  are  zero.  These  parameters 
are  somewhat,  but  not  completely,  unlike  those  of  an  F- 
18  trimmed  at  the  corresponding  flight  conditions. 


Table  1:  Effector  Rate  and  Position  Limits 


Effector 

Position  Limits 
(deg) 

Rate  Limits 
(deg) 

Horizontal 

Tails 

+10.5,  -24 

±40 

Ailerons 

+45,  -25 

±100 

Leading 
Edge  Flaps 

+33,  -3 

±15 

Trailing 
Edge  Flaps 

+45,  -8 

±18 

Rudder 

+30,  -30 

±82 

Control  Law  Designs 

For  comparative  purposes,  two  adaptive  control  law 
designs  and  one  non-adaptive  design  will  be 
implemented.  The  first  adaptive  design  includes 
estimation  of  the  full  set  of  17  control  effectiveness 
parameters.  The  second  adaptive  design  entails  defining 
a  set  of  three  pseudo  controls  and  then  estimating  only 
five  pseudo  control  effectiveness  parameters. 

Full  Control  Estimation 

This  design  essentially  follows  that  of  Reference  8  and 
is  similar  to  that  of  References  7  and  9.  However,  the 


2 

American  Institute  of  Aeronautics  and  Astronautics 


Table  2;  Aircraft  Model  Parameters  for  trimmed  flight  at  30,000  ft,  Mach  0.7  (all  angles  and  angular  rates 
are  assumed  to  be  in  radians  and  radians/sec  respectively) 


=-11.04 

II 

O 

/,  =0.4164 

//to  =-19.72 

/,„  =4.709 
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V 
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II 
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n.  =-0.0698 

rig^  =  -0.0698 

n.  ^-\J422 

F= 696  ft/sec 

«(,  =0,0681 

6»o  =0.0681 

Table  3;  Aircraft  Model  Parameters  for  trimmed  flight  at  40,000  ft,  Mach  0.6  (all  angles  and  angular  rates  are  assumed 
to  be  in  radians  and  radians/sec  respectively) 
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following  design  does  not  contain  the  fuzzy  logic 
mitigation  of  actuator  saturation  as  was  done  in 
Reference  8.  The  aircraft  equations  of  motion  (Eq.  1) 
can  be  put  in  the  following  form: 

j  =  Oo  (x, )  +  <D(x,  )w,  +  5(x,  )ty 
^  =  (^)  +  ¥\  +  D{x,  )u  (2) 

»7  =  9o  W  +  (x)w  +  (x,  )u 

In  this  form,  is  a  vector  of  the  outputs  that  will  be 
controlled,^  is  the  virtual  control  vector  used  to 
control jp,  7]  is  the  vector  of  uncontrolled  states, 
>Vj,W2,  and  are  the  vectors  of  unknown 

parameters,  is  a  subset  of  the  state  vector  x ,  and  U 
is  the  vector  of  control  effector  commands  where, 

y  =  {(j),a,pY 

<i)^{p,q,rY  (3) 

r)  =  6 

X,  =  {a,p,^,9) 


where  Kq  is  a  positive  definite  symmetric  matrix  and 
is  the  output  of  a  command  generator  that  is  a 
linear,  stable,  3rd  order  system.  It  will  be  useful  to 
define  the  derivative  of  S  as: 

s  =  <I>(Xi  )Wj  +  B(x^  )(a)j  +  5)  +  V 
co-a)j  -¥6) 
v  =  <t>o-y,+K^e 

In  Equation  5,  0)^  represents  the  desired  value  of  the 
virtual  controls.  However,  since  these  are  states  and 
not  directly  controllable  effectors,  there  will  be  an  error 
0). 

We  will  next  choose  0)^  through  the  use  of  the 
following  Lyapunov  function: 

U,  =  (s^s +  L^w^)/2 

where  ij  is  a  positive  definite  diagonal  matrix  and 
is  the  estimate  of  . 

Taking  the  derivative  of  t/j  yields: 

U^  Z,Mi  (7) 


We  will  next  define  the  error,  e  ,and  a  function  s  that 
combines  error  and  integrated  error, 
e  =  y-y, 

s  =  e  +  KqX^  (4) 

Xs  =e 


0)^  can  then  be  chosen  as: 

=  B-'  (x^  ,5  -  (x,  )m',  -  v(x,  )]  (8) 


where  is  a  positive  definite  diagonal  matrix. 
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As  a  result  of  that  choice  of  Q)^ , 

f/j  + +  (9) 

If  we  put  <2?^  in  the  following  form: 

(i><i=¥od-¥xd^y-^¥id^2  00) 

then, 


Following  the  approach  of  Reference  7,  it  can  be  shown 
that  for  the  nominal  system,  S  ,  and  e  all  tend  to 
zero  as  f  — >  00 . 

Furthermore,  it  is  desirable  to  bound  the  parameter 
estimates  on  some  interval  to  help  avoid  singularities  in 
the  control  law.  To  this  end,  a  projection  algorithm  is 
employed  similar  to  that  described  in  Reference  9. 
Suppose  that  each  component  of  (i  =  1,  2,  u)  is 
known  to  lie  in  an  open  interval  Wik  Then  a 

modified  adaptation  law  is  chosen  as  (i  =  1 ,  2,  u), 


^  =  ¥oa+¥u^\-^¥A-^  D(x,  w„  )u 

¥oa  =¥o -¥od  +¥20^1 

¥xa=¥u 

¥2a=¥x-¥ld 

where  =  w,.  from  equation  15. 

We  will  next  chose  a  2nd  Lyapunov  function  such  that, 


if  Wf*  e  W A  or  w,*  =  m,*  and  <  0 
orM-,*=4and/i,,>0 

0  otherwise 


(18) 


IJ.^=IJ^  +(S^5  + +wILjW^)I2  (12) 

where  L^  and  are  positive  definite  diagonal 
matrices.  Taking  the  derivative  of  gives: 

=-//i:,,s+wf[<D^s+I,w,  wJLjWj  (-13) 

+Dt\ 

We  will  next  chose  the  following  control  law: 

bu  =  £)(x,  w„  )u  =  (-B^s  -  y/oa  “  ^22^)  =  '”*1  (i^) 

and  the  following  parameter  update  laws: 

w,  =Z7'(<D^5  +  ^f» 

^2=^'¥L^ 

where  is  chosen  such  that, 

[D{x,  )  -  D(x,  )]u  =  'V„{x,u)w„  (1 6) 

As  a  result,  since  AT,,  and  K22  are  positive  definite 
diagonal  matrices, 

U2=-s^  K22a<0  (17) 


Pseudo  Control  Estimation 

A  separate  adaptive  control  law  can  be  obtained  if  we 
define  the  following  pseudo  controls  (Up)  and 
corresponding  control  effectiveness  (Dp): 


0 

Ig  '] 

Odir 

(8  \ 

D.u„  = 

p  p 

0 

nix 

0 

^lon 

with. 

0 

=(4  Jji  Y  (20) 

The  control  law  design  then  proceeds  exactly  as 
described  for  Full  Control  Estimation  except  that  Dp,  Up, 
and  Wup  should  replace  D,  w,  and  Wu  in  equations  2 
through  18.  In  order  to  determine  the  actual  control 
commands  from  the  pseudo  control  commands,  another 
relationship  is  required.  Namely, 

£>m  =  (21) 

where  D  now  indicates  an  estimate  of  the  full  control 
effectiveness  matrix  that  is  either  given  a  priori  or 
generated  on-line  from  a  source  other  than  the  nonlinear 
adaptive  controller. 

Non- Adaptive 

A  non-adaptive  controller  that  is  essentially  a  two-loop 
dynamic  inversion  controller  is  obtained  by  simply 
turning  off  the  parameter  updates  ( w,.  =  0  (/  =  1,2,  u)) 
for  either  of  the  preceding  adaptive  designs. 
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Control  Allocation  Techniques 


Discrete  Time  Direct  Allocation 


In  the  above  control  law  derivations,  the  control  vector 
(m)  is  chosen  such  that  either, 


Ax7W7x1  =  f»dcs, 


(22) 


or, 


Ax7"7x1  -^p,a 


3x3 


(23) 


where  the  subscripts  (m  x  n)  represent  the  matrix  and 
vector  dimensions  for  the  chosen  design  of  Equations  1 
and  2.  For  Full  Control  Estimation  (Equation  22),  the 
elements  of  Dare  updated  by  the  parameter  update  law 
of  Equation  15.  For  Pseudo  Control  Estimation 
(Equation  23),  it  is  the  elements  of  D^that  are  updated 

by  the  parameter  update  law  of  Equation  15.  When 
there  are  no  constraints  on  the  control  vector,  Equations 
22  and  23  will  generally  have  an  infinite  number  of 
solutions.  However,  when  control  constraints  are 
considered,  these  equations  can  either  have  an  infinite 
number  of  solutions,  no  solution,  or  a  unique  solution. 
The  solution  of  Equation  22  or  23  for  the  control  vector 
(m)  is  the  control  allocation  problem.  Only  a  brief 
description  of  each  control  allocation  algorithm  will  be 
presented.  For  more  detailed  information,  consult  the 
references  cited. 


Direct  Allocation 


Durham’®  recognized  that  existing  control  allocation 
schemes  did  not  fully  exploit  the  capabilities  of  the 
redundant  control  effectors  they  were  commanding. 
Durham  therefore  defined  the  Attainable  Moment 
Subset  (AMS)  to  be  the  set  of  all  moments  capable  of 
being  produced  by  a  given  set  of  control  effectors 
constrained  to  move  only  within  certain  bounds.  In 
Direct  Allocation  (DA),  the  control  effector  commands 
are  calculated  as  a  scaled  version  of  the  deflections  that 
yield  the  maximum  attainable  moment  in  the  desired 
direction.  That  is,  given  the  direction  of  any  rrides,  the 
controls  are  found  that  lie  on  the  boundary  of  the  AMS 
in  that  direction.  If  the  desired  moment  is  smaller  than 
maximum  attainable  moment  in  the  given  direction,  the 
control  vector  is  simply  scaled  down  (preserving  the 
direction)  such  that  Equation  22  or  23  is  satisfied.  If  the 
desired  moment  is  greater  than  or  equal  to  the 
maximum  attainable  moment,  then  the  control  vector  is 
set  equal  to  that  on  the  boundary  of  the  AMS.  Direct 
allocation  insures  that  all  moments  attainable  by  the 
control  effectors  are  indeed  realizable  by  the  control 
allocation  algorithm.  For  those  moments  outside  the 
AMS,  the  direction  of  the  moment  is  preserved. 


Direct  Allocation  only  considers  position  limiting  of  the 
control  effectors.  A  straightforward  discrete  time 
extension  of  Direct  Allocation  allows  the  incorporation 
of  effector  rate  limiting  to  be  considered  as  well.  In 
Discrete  Time  Direct  Allocation”'  (also  referred  to  as 
Moment-Rate  Allocation)  the  problem  essentially  can 
be  defined  as  the  solution  to, 


=  (24) 

Now,  the  limits  on  the  control,  Aw,  are  taken  as  the 
more  restrictive  of  1)  the  global  position  limits  or  2)  the 
effector  rate  limits  multiplied  by  A/.  That  is,  the  control 
effectors  are  limited  by  how  far  they  can  travel  during 
the  next  time  increment  (Ar).  Therefore,  the  AAMS  can 
be  defined  as  the  set  of  all  moments  producible  over  the 
next  time  step.  With  this  problem  formulation,  the  same 
algorithms  used  in  Direct  Allocation  can  be  applied  to 
Equation  24. 

Discrete  Time  Direct  Allocation  (DTDA)  is  prone  to  a 
form  of  controller  windup”'  To  alleviate  this 
problem  various  restoring  techniques  have  been  offered 
which  have  the  ability  to  drive  the  control  allocation 
solution  so  as  to  optimize  a  given  objective  function. 
The  results  presented  herein  use  the  minimum  norm 
restoring  technique  as  described  in  References  1 1  and 
12.  This  technique  drives  the  solution  to  that  which 
minimizes  the  2-norm  of  the  control  vector. 

Pseudo-Inverse  and  Weighted  Pseudo-Inverse 

A  simple  and  computationally  efficient  solution  to  the 
control  allocation  problem  is  the  use  of  a  Weighted 
Pseudo-Inverse  (WPI).  In  particular,  the  solution  of 
Equation  22  or  23  that  minimizes  the  following 
quadratic  function: 


J  =  u^Wu  (25) 

where  W  is  an  appropriate  weighting  matrix,  is  given 
by‘^ 

u  =  W-\by[bW-\bYrm^,  (26) 

While  Equation  26  provides  a  unique  solution  to 
Equation  22  or  23,  it  does  not  consider  any  limitations 
on  the  control  effectors.  This  solution  technique  is  used 
as  the  sole  control  allocator  in  a  similarly  designed 
adaptive  controller  in  Reference  9.  When  W=I,  the 
solution  is  termed  the  Pseudo-Inverse  (PI)  solution. 
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Simulation  Results 

The  aircraft  simulation  model  is  defined  in  Equation  1. 
The  control  effector  rate  and  position  limits  are  given  in 
Table  L  The  aircraft  model  parameters  for  two  separate 
flight  conditions  are  given  in  Tables  2  and  3.  The  three 
control  laws  are  defined  as  specified  previously  with 
the  controller  gains  as  follows: 

Ko=0.4I,  K, i-I,  K22-4I,  L,-0.2I,  L2=0.2I,  L3=0.2I  (27) 

where  I  is  the  identity  matrix  of  the  proper  dimension. 

The  initial  condition  estimates  of  the  stability  and 
control  parameters  for  the  Full  Control  Estimation 
design  were  chosen  to  vary  ±50%  from  their  nominal 
values  as  indicated  in  Table  4.  For  the  Pseudo  Control 
Estimation  design,  the  initial  pseudo  control 
effectiveness  matrix  is  chosen  as  the  identity 


Table  4;  Initial  Conditions  for  Parameter  Estimates 


IC  Number 

w,  (i  =  l,2) 

1 

O.SOWj 

0.50w„ 

2 

1  .OOw, 

0.50^^ 

3 

1.50w, 

0.50w„ 

4 

0.50w, 

1.00>v„ 

5 

1  .OOw, 

l-OOw^ 

6 

1.50w^ 

1.00w„ 

7 

O.SOh', 

1.50w„ 

8 

I.OOm', 

1.50m^„ 

9 

1.50w, 

matrix  and  the  estimated  control  effectiveness  matrix 
(D)  is  fixed  at  the  value  determined  by  the  control 
effectiveness  parameter  estimates  given  in  Table  4.  The 

parameter  estimates,  (i  =  1,  2,  u),  are  bound  by  the 
projection  algorithm  of  Equation  18  to  be  within  75% 
of  their  nominal  values  with  the  exception  of 

L  and/zc  which  are  bound  between -0.2  and  0.2. 

Both  medium  and  large  amplitude  single  and  multi-axis 
doublet  maneuvers  were  simulated  at  both  flight 
conditions  starting  with  the  initial  estimates  given  Table 
4.  The  piecewise  constant  inputs  to  the  command 
generator  were  of  the  form:  y*(t)  =  ((j)/,  ao+ai*,  p/)  for 
0  <  t  <  4  seconds,  y*(t)  =  ((t)2*,  ao+a2*,  P2*)  for  4  <  t  <  8 
seconds,  and  y*(t)  =  (0,  ao,  0)  for  8  <  t  <  15  seconds. 
Table  5  lists  the  command  generator  inputs  for  the  six 


maneuvers  considered.  Unless  otherwise  noted,  the 
initial  conditions  chosen  for  the  command  generator 
=  >,(0)  =  y^(0)  =  0  ,  except  ac(0)=ao,  and 
the  command  generator  poles  are  placed  at  -1.5  and 
-1.5±1.06i. 


Table  5:  Command  Generator  Inputs 


Maneuver 

Number 

^t'lXdeg) 

<l>2’(deg) 

ar(deg) 

aa’Cdeg) 

1 

45 

-45 

0 

0 

2 

90 

-90 

0 

0 

3 

0 

0 

20 

-10 

4 

0 

0 

40 

-20 

5 

30 

-30 

15 

-7.5 

6 

60 

-60 

30 

-15 

The  only  difference  among  the  simulations  for  each 
control  allocation  method  is  the  algorithm  used  to  solve 
Equation  22  or  23  for  the  control  vector.  The  algorithms 
used  for  both  Direct  Allocation  and  Discrete  Time 
Direct  Allocation  were  adapted  from  Reference  14. 
Weighted  Pseudo-Inverse  allocation  is  given  explicitly 
by  Equation  26  where  the  weighting  matrix  is  chosen 
as,  W=  diag(30,  30,  3,  3,  8,  9,  4).  Simulation  results  are 
also  presented  for  the  Pseudo  Inverse  (PI)  solution 
(i.e.,W=I). 

Simulations  of  all  six  maneuvers  were  carried  out  at 
both  flight  conditions  for  all  nine  initial  conditions, 
three  control  law  designs,  and  with  each  of  the  four 
control  allocation  techniques.  All  simulations  were  run 
for  15  seconds.  Since  the  maneuvers  of  Table  5  should 
be  completed  within  12  seconds,  it  is  of  interest  to  note 
how  many  individual  simulation  runs  ended  with  large 
tracking  errors  (average  norm  of  the  tracking  error 
greater  than  10  degrees)  between  13  and  15  seconds. 
The  existence  of  large  tracking  errors  during  the  last 
two  seconds  is  used  as  an  indication  the  controller  was 
unable  to  perform  adequate  tracking  for  the  specified 
maneuver. 

Table  6  shows  the  number  of  simulations  terminating 
with  large  tracking  errors  for  each  control  allocation 
technique.  If  the  data  of  Table  6  is  broken  down  further, 
it  is  found  that  the  runs  ending  with  large  tracking  error 
occur  for  two  distinct  situations.  For  the  non-adaptive 
controller,  the  entries  in  Table  6  all  occur  at  the  30,000 
ft.,  Mach  0.7  flight  condition.  On  the  other  hand,  the 
entries  in  Table  6  for  both  adaptive  designs  all  occur  at 
the  40,000  ft.,  Mach  0.6  flight  condition. 

Of  the  cases  indicated  Table  6  for  the  non-adaptive 
design,  the  vast  majority  occur  for  initial  conditions 
where  the  control  effectiveness  is  underestimated. 
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Table  6:  Number  of  Cases  Terminating  with  Large 
Tracking  Errors _ 


DA 

PI 

WPI 

DTDA 

Non-Adaptive 

5 

5 

6 

5 

Full  Control 
Estimation 

2 

1 

0 

0 

Pseudo  Control 
Estimation 

9 

5 

0 

9 

Therefore,  it  is  likely  that  the  control  law/control 
allocator  combination  over  commands  the  actuators 
leading  to  poor  damping  and  large  tracking  errors.  The 
data  in  Table  6  reveals  that  only  three  cases  yield  large 
terminal  tracking  errors  for  the  Full  Control  Estimation 
design  (2  for  Direct  Allocation  and  1  for  Pseudo- 
Inverse  allocation).  The  Pseudo  Control  Estimation 
design,  however,  has  significant  trouble  at  the  lower 
dynamic  pressure  flight  condition.  A  notable  exception 
is  the  Weighted  Pseudo-Inverse  allocation  method  that 
had  no  cases  ending  with  large  tracking  errors  for  either 
of  the  adaptive  designs. 

In  order  to  present  meaningful  tracking  error  statistics, 
two  reduced  data  sets  will  be  used.  The  first  data  set 
results  from  removing  the  nine  distinct  cases  that  lead 
at  least  one  allocation  scheme  to  terminate  with  large 
tracking  errors  for  either  the  non-adaptive  design  or  the 
Full  Control  Estimation  design.  The  word  “case”  in  this 
context  refers  to  a  particular  combination  of  flight 
condition,  maneuver  number,  and  initial  condition 
number.  The  results  for  the  average  2-norm  of  the 
tracking  error  over  all  maneuvers  in  this  data  set  are 
given  in  Table  7.  The  data  in  Table  7  shows  that  all 
allocation  techniques  perform  reasonably  well  for  the 
non-adaptive  design  with  pseudo-inverse  allocation 
winning  by  a  narrow  margin.  However,  when  Full 
Control  Estimation  adaptation  is  employed,  Weighted 
Pseudo- Inverse  allocation  vaults  into  the  lead  followed 
closely  by  Discrete  Time  Direct  Allocation.  The  values 
for  the  Pseudo  Control  Estimation  design  are  presented 
in  Table  7  to  illustrate  the  average  norm  of  the  tracking 
error  is  reduced  by  30%  in  going  from  the  Full  Control 
Estimation  to  the  Pseudo  Control  Estimation  control 
law  design  for  Weighted  Pseudo-Inverse  allocation. 
This  represents  a  reduction  of  over  60%  from  the 
smallest  average  tracking  error  norm  for  the  non- 
adaptive  design.  The  Pseudo  Control  Estimation  results 


given  in  Table  7  for  the  remaining  control  allocation 
techniques  contain  data  from  maneuvers  ending  with 
large  tracking  errors  and  are  hence  not  directly 
comparable. 

The  second  data  set  results  from  removing  the  eleven 
additional  cases  that  also  lead  at  least  one  allocation 
scheme  to  end  with  large  tracking  errors  for  the  Pseudo 
Control  Estimation  design.  The  average  norms  of  the 
tracking  errors  for  this  data  set  are  given  in  Table  8. 
The  data  in  Table  8  illustrates  the  same  trends  as  the 
data  in  Table  7  with  the  added  information  that,  when  it 
works,  Pseudo  Control  Estimation  outperforms  Full 
Control  Estimation  for  each  of  the  control  allocation 
technique  tested. 

On  the  surface,  these  results  are  quite  surprising  since 
the  given  Weighted  Pseudo-Inverse  allocation 
technique  produces  admissible  solutions  for  only  about 
10%  by  volume  of  the  AMS.  An  admissible  solution  is 
defined  as  one  for  which  the  estimated  attained  control 

effect  )  equals  the  desired  control  effect  (mjes  or 
Up).  Here,  Uum  represents  the  controls  as  constrained  by 
their  position  limits.  Perhaps  more  importantly, 
weighted  pseudo-inverse  allocation  also  produces 
admissible  solutions  for  only  10%  by  volume  of  the 
AAMS  as  defined  by  the  control  effector  rate  limits.  By 
contrast  (and  by  definition)  Discrete  Time  Control 
Allocation  produces  admissible  solutions  for  100%  of 
the  AMS  as  well  as  for  100%  of  the  AAMS, 

A  deeper  understanding  of  the  contributing  control 
law/control  allocation  interactions  can  be  gained  by 
detailed  examination  of  a  single  maneuver.  For  this 
examination,  a  360°  rolling  maneuver  at  constant  angle- 
of-attack  and  zero  sideslip  will  be  commanded  at  the 
lower  dynamic  pressure  flight  condition  of  Table  3.  The 
inputs  to  the  command  generator  are  y*(t)  =  (90°t,  Wq, 
0)  for  0  <  t  <  4  seconds  and  y*(t)  (360°,  ao,  0)  for  4  < 

t  <  15  seconds.  The  command  generator  poles  for  this 
maneuver  are  placed  at  ^.5  and  ~4.5±3. 1 8i. 

Figures  la-lh  show  the  simulation  results  for  the 
Pseudo  Control  Estimation  control  law  with  WPI 
control  allocation  where  the  initial  parameter  estimates 


Table  8:  Average  Error  Norm  for  Reduced 
Maneuver  Set  (Twenty  Cases  Removed) 


DA 

PI 

WPI 

DTDA 

Non-Adaptive 

2.6145 

2.3670 

2.5135 

2.5118 

Full  Control 
Estimation 

1.4043 

1.4567 

1.1689 

1.2404 

Pseudo  Control 
Estimation 

0.7439 

1.0214 

0.7257 

0.8352 

Includes  cases  with  large  tracking  errors  at  termination 
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are  taken  to  correspond  to  IC  Number  1  from  Table  4. 
In  Figures  la  and  lb,  the  simulated  roll  angle  and 
angle-of-attack  (solid  lines)  are  seen  to  track  the 
commanded  values  (dotted  lines)  acceptably  well.  The 
tracking  errors  are  explicitly  plotted  in  Figure  Ic. 
Skipping  to  Figures  le  and  If,  the  control  surface 
deflections  exhibit  little  rate  or  position  saturation  with 
the  exception  of  the  rudder.  The  rudder  is  either  rate  or 
position  limited  during  most  of  the  maneuver. 

The  changes  in  the  parameter  estimates  during  the 
maneuver  are  shown  in  Figure  Ig.  Though  these 
changes  are  mostly  smooth,  the  effects  of  the  parameter 
bounding  algorithm  can  be  seen  to  kick  in  during 
periods  of  rudder  position  saturation.  Figure  Id  shows 
the  pseudo  control  command  and  Figure  Ih  gives  the 
difference  between  the  desired  pseudo  control  and  the 
estimated  control  effect  produced  by  the  control 
allocator.  The  nonzero  values  of  Figure  Ih  indicate  that 
the  desired  pseudo  control  is  outside  the  AAMS  for  the 
control  allocator. 

Figures  2a-2h  show  simulation  results  for  the  same 
situation  as  Figure  1,  except  that  DTDA  is  used  as  the 
control  allocation  technique.  Figures  2a  -  2c  show  the 
tracking  performance  to  be  unacceptable.  Figures  2e 
and  2f  show  severe  limiting  of  control  surfaces.  The 
behavior  of  the  control  law  during  periods  of  control 
surface  position  limiting  is  one  of  the  key  interactions 
responsible  for  the  relatively  poor  results  for  DTDA. 
Consider  the  system  behavior  between  2  and  4  seconds. 
During  this  time,  five  control  surfaces  reach  their 
position  limits  indicating  the  desired  control  effect  lies 
either  outside  the  AMS  or  on  the  boundary  of  the  AMS. 
With  the  controls  this  severely  limited,  the  system  does 
adequately  respond  to  changes  in  the  desired  control 
effect  and  hence  the  parameter  estimates  all  head  to 
their  limiting  values.  It  is  exactly  this  type  of  behavior 
that  lead  to  the  use  of  fuzzy  logic  to  reduce  the 
adaptation  and  controller  gains  when  the  controls 
become  limited^.  As  the  parameter  estimates  become 


saturated,  the  controller  looses  damping  and  the 
response  is  greatly  degraded. 

To  understand  the  other  contributing  interaction,  it  is 
necessary  to  examine  the  difference  between  the  pseudo 
control  commands  and  the  estimated  control  effect. 
Table  9  shows  this  information  for  the  first  0.04 
seconds  of  the  maneuver  for  both  WPI  and  DTDA 
control  allocation.  For  the  first  0.02  seconds,  both 
allocation  techniques  produce  admissible  controls.  At 
0.03  seconds,  the  DTDA  allocator  still  produces 
admissible  controls  but  the  WPI  allocator  does  not. 
Finally,  neither  allocator  produces  admissible  solutions 
at  0.04  seconds. 

The  important  interaction  lies  with  the  how  the  control 
allocation  routines  deal  with  unattainable  values  of  the 
pseudo  control  commands.  While,  DTDA  preserves  the 
direction  of  the  desired  effect  by  design,  WPI  can 
produce  unexpected  results.  At  0.04  seconds,  the 
DTDA  routine  provides  87%  of  the  desired  control 
effect  while  the  WPI  solution  provides  104%  of  the  5iat 
desired  and  only  73%  of  the  desired  5dir.  The  WPI 
solution  sacrifices  directional  control  power  and 
actually  yields  more  lateral  control  than  desired.  Since 
the  WPI  solution  generates  artificial  lead  in  the  lateral 
axis,  the  magnitudes  of  the  pseudo  control  commands 
are  actually  reduced  while  performing  the  given  rolling 
maneuver  and  the  stability  of  the  control  law  is 
enhanced. 

However,  the  lead  produced  by  the  WPI  solution  could 
serve  to  destabilize  the  control  law  for  another 
maneuver.  For  example,  if  the  command  generator 
inputs  are  taken  to  be  y*(t)  =  (35°sin(2t),  Oq,  0)  for  0  <  t 
<  10,  then  the  roll  angle  responses  shown  in  Figure  3 
result.  This  time,  the  DTDA  technique  exhibits 
adequate  tracking  while  the  WPI  solution  generates 
large  errors. 


Table  9:  Difference  Between  Pseudo  Control  commands  and  Estimated  Control  Effect 


Desired  Control  Effect  (Up) 

Desired  Control  Effect  -  Estimated 

Control  Effect  for  WPI  (Up-  ) 

Desired  Control  Effect  -  Estimated 

Control  Effect  for  DTDA  (Up-  ) 

Time 

(sec) 

■ 

m 

0.01 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

ISSS2USI 

0.0000 

0.0000 

itMii 

KtliTdil 

0.0000 

0.0000 

0.04 

■lima 

0.0000 

0.0171 

0.0000 

0.0000 

0.0022  1 

8 

American  Institute  of  Aeronautics  and  Astronautics 


Time  (sec)  Time  (sec) 


Figure  1 :  Simulation  Results  for  WPI  Allocation 
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a  (deg) 


0  2  4  6  8  10  0  2  4  (h)  6 

Time  (sec)  Time  (sec) 

Figure  2:  Simulation  Results  for  DTDA  Allocation 
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Figure  3:  Roll  Angle  Responses  for  Sinusoidal  Input 


Conclusions 

The  effect  of  the  choice  of  control  allocation  algorithm 
used  in  conjunction  with  a  nonlinear  adaptive  control 
law  has  been  studied.  The  results  show  that  this  choice 
can  have  dramatic  impact  on  system  performance.  For 
the  set  of  maneuvers  and  initial  conditions  tested,  the 
Weighted  Pseudo-Inverse  solution  performed  well 
despite  producing  admissible  solutions  for  only  10%  of 
the  Attainable  Moment  Subset. 

The  less  admirable  system  performance  that  results 
from  using  Disrcrete  Time  Direct  Allocation  with  the 
backstepping  adaptive  designs  is  attributed  to  an 
important  control  law/control  allocation  interaction  that 
results  when  the  effectors  become  rate  or  position 
saturated.  Since  the  situation  could  have  been  quite 
different  if  the  actuator  mitigation  logic  of  Reference  8 
had  been  included,  these  results  must  not  be  viewed  as 
a  detriment  to  Direct  Allocation.  However,  the  results 
do  serve  to  illustrate  that  simply  using  the  best  (in  terms 
of  ability  to  meet  attainable  moment  demands)  control 
allocator  does  not  necessarily  yield  the  best  system 
performance.  To  further  complicate  matters,  the 
performance  results  can  be  sensitive  to  changes  in  the 
maneuver  set,  the  flight  conditions,  and  the  initial 
parameter  estimates.  Hence,  the  results  presented  herein 
are  meant  to  illustrate  the  important  control  law/control 
allocation  interactions  rather  than  to  judge  the  relative 
merits  of  the  control  laws  or  control  allocation 
techniques  on  an  individual  basis. 

A  better  understanding  of  some  of  the  basic  control 
law/control  allocation  interactions  facilitates  an 
integrated  design  that  endeavors  to  attain  the  best 
performance.  Integration  of  an  extended  set  of  control 


allocation  schemes  with  a  broad  array  of  nonlinear 
control  designs  is  planned  for  future  research. 
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